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ABSTRACT 



We present an analysis of elemental abundances of He, N, O, Ne, S, and Ar in 

Ph Magellanic Cloud planetary nebulae (PNe), and focus initially on 14 PNe in the Small 

Magellanic Cloud (SMC). We derived the abundances from a combination of deep, high 

dispersion optical spectra, as well as mid-infrared (IR) spectra from the Spitzer Space 

Telescope. A detailed comparison with prior SMC PN studies shows that significant 



o 

?-H variations among authors of relative emission line flux determinations lead to systematic 

^ discrepancies in derived elemental abundances between studies that are > 0.15 dex, in 

'— ' spite of similar analysis methods. We used ionic abundances derived from IR emission 

^_i lines, including those from ionization stages not observable in the optical, to examine 

^ the accuracy of some commonly used recipes for ionization correction factors (ICFs). 

pg These ICFs, which were developed for ions observed in the optical and ultraviolet, 

^N relate ionic abundances to total elemental abundances. We find that most of these 

ICFs work very well even in the limit of substantially sub-Solar metallicities, except for 

(^ PNe with very high ionization. Our abundance analysis shows enhancements of He and 

^^ N that are predicted from prior dredge- up processes of the progenitors on the AGB, as 
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well as the well known correlations among O, Ne, S, and Ar that are little affected by 
nucleosynthesis in this mass range. We identified MG 8 as an interesting limiting case 
of a PN central star with a ~ 3.5 Mq progenitor in which hot-bottom burning did not 
occur in its prior AGB evolution. We find no evidence for O depletion in the progenitor 
AGB stars via the 0-N cycle, which is consistent with predictions for lower-mass stars. 
We also find low S/0 ratios relative to SMC H II regions, with a deficit comparable to 
what has been found for Galactic PNe. Finally, the elemental abundances of one object, 
SMP-SMC 11, are more typical of SMC H II regions, which raises some doubt about its 
classification as a PN. 

Subject headings: Magellanic Clouds — planetary nebulae: general - stars: evolution 



Introduction 



Studies of planetary nebulae (PNe) in the Magellanic Clouds over the past two decades have 
lead to substantial progress in understanding the late stages of stellar evolution, particularly the 
nucleosynthesis and mixing processes in their progenitor AGB stars that lead to the enrichment of 
He, C, and N in the interstellar medium. PNe in the Magellanic Clouds (MCPNe) are well suited to 
this endeavor in that their average metallicity, being substantially lower than Solar, constrains and 
challenges stellar evolution theory in ways not easy to test in other environments. Moreover, the 
progenitors of MCPNe span a broad range of abundances, owing to both the large range in main- 
sequence masses (roughly 1-8 Mq , with a correspondingly large range of pre- AGB ages) , and to the 
increase in the average abundances in the ISM of the host galaxy over time, mostly from the ejecta 
of supernovae. MCPNe are easy to identify because of their bright nebular emission lines, they are 
numerous, and on average their emission is little affected by interstellar extinction. With the aid 
of high-resolution imaging, MCPNe are usually resolved and both the nebula and its central star 
can be easily observed against an often complex background of stellar and diffuse nebular emission. 
This set of observables and a known distance is a combination that cannot be matched for any other 
sample with today's technology. A good deal of spectroscopic data have been published on MCPNe 
in the past two decades, mostly in the optical, by Monk, Barlow & Clegg' (1988), Meatheringham 
& Dopita ( 1991a|[b ), Vassiliadis, et al. (1992), and Leisy fc Dennefeld| ( |1996 who also included 



UV spectroscopy); these data were compiled and augmented by Leisy & Dennefeld (2006). More 



recently Stanghellini, et al. (2002), Stanghellini, et al. (2003), and Shaw, et al. (2006) published 



HST /STIS slitless optical and UV spectra; additional optical spectroscopy was published by Idiart, 



Maciel & Costa (2007); and Stanghellini, et al. (2009) published carbon abundances based on UV 



spectroscopy. In all, abundances in more than 180 MCPNe have been derived and analyzed. 

The most comprehensive analysis of abundance patterns among MCPNe to date is that of 



Leisy & Dennefeld (2006 hereafter, LD06), who explored the distinction between PNe of Type I 
vs. non-Type I (Peimbert||1978 Torres-Peimbert fc: Peimbert|199T), the possibility of O destruction 



by nucleosynthesis in the core of the progenitor AGB star via the 0-N cycle, the efficiency of the 
third dredge-up at low metallicity, and the constancy of Ne, Ar, and S through AGB evolution, 



among other relationships. More recently, Maciel, Costa & Idiart (2009) also examined abundance 



patterns in MCPNe, using many of the above-cited data in addition to new observations by Idiart, 



Maciel &: Costa (2007, hereafter, IMC07). The various reported trends in elemental abundances in 
both studies had a good deal of scatter, which complicated the comparisons with predictions from 
stellar evolution theory for elemental yields. Indeed, LD06 noted that the spectrophotometry taken 
from the literature for several PNe was of rather poor quality, either because critical faint emission 
lines were not detected, or because the discrepancies with their own data were very large. In the end 



LD06 excluded many of the emission line fluxes from Monk, Barlow & Clegg (1988) and Vassiliadis, 



et al. ( 1992 ) as unreliable; for the remaining PNe they attempted to remove some common sources 



of discrepancy between the various published abundances by re-deriving the abundances using their 
own determination of the nebular electron temperature (Tg) and density {N^), a consistent set of 
supporting atomic data, and using a uniform set of prescriptions for ionization correction factors 
(ICFs). LD06 estimated average uncertainties in their final elemental abundances in the range 
0.05 dex for He to as high as 0.5 dex for S. 

High accuracy of elemental abundance determinations from PNe is crucial for validating stellar 
evolution models, particularly the chemical yields, to useful precision. Achieving high accuracy de- 
mands great care in the observations and their calibration, a robust derivation of the abundances, 
and a thorough accounting of the observational uncertainties. Yet the subject of observational 



uncertainties is seldom discussed in quantitative depth in the literature, although Stasihska (2004) 



provides a fairly complete qualitative description of many of the key issues. A comparison of abun- 
dances of elements in individual nebulae among the above-cited papers shows differences as large 
as 1.0 dex or more, with differences of ~ 0.15 dex being very common. This level of disagreement 



among authors generally has been noted before (for MCPNe specifically by, e.g., LD06, and Maciel, 
Costa &: Idiart||2009 ), and in some cases it is recognized as a significant problem (see, e.g., Pottasch 



& Beintema 1999) for comparing observations to chemical yield predictions from theory. Close 



examination of the differences reveals a variety of causes, including discrepant emission lines fluxes 
(because of calibration errors, limited signal-to-noise, blending with other, closely spaced emission 
features, etc.), differences in adopted extinction constants, different adopted physical parameters 
(particularly Tg), different adopted atomic data (particularly collision strengths), and the adoption 
of different ICF formulae. These systematic discrepancies can form a cascade, where errors at each 
stage of the abundance calculation can magnify errors in subsequent stages in the chain. 

We undertook a program of spectrophotometry of MCPNe with the aim of improving the 
accuracy of abundance determinations to support a more precise comparison of the chemical yields 
that are predicted by stellar evolution theory. We are able to improve upon prior results by 
leveraging highly accurate, space-based emission line fluxes of brighter emission lines, adding new 
mid-IR fine structure emission lines from some ions (which are insensitive to Tg, including those 
of [O IV], [Ne II], [Ne III], [Ne V], [S III], and [S IV]), and by obtaining additional, high-resolution 



ground-based optical spectra. In the process we have characterized many systematic difficulties of 
prior studies, and we have endeavored to quantify our uncertainties with great care. We focus in 
this paper on 14 PNe in the Small Magellanic Cloud (SMC). We describe the observing programs, 
data reduction, and emission line measurements in §2. In §3 we describe the derivation of the 
nebular physical diagnostics, and of ionic and elemental abundances, and compare them with those 
in the literature in detail. In §4 we probe the systematic discrepancies and possible sources of error 
in elemental abundances between authors, and analyze our adopted abundances in the context of 
stellar evolution in a low metallicity environment. We conclude in §5 with a summary and a few 
thoughts on future studies of abundances in emission line nebulae. 



2. Observations 

Our strategy for optical observations was to obtain spectra of our target PNe that were deeper, 
had higher resolution, and broader wavelength coverage than were available in the literature. There 
are several advantages to this approach: Higher spectral resolution allows us to resolve closely 
spaced lines such as the [O II] pair at A3726.0, 3728.8; [Ne III] A3967.4 and He A3970.1; [Ar IV] 
A4711.3, He I A4713.1 and Ne IV A4715.4; He II A4859.3 and H/3 A4861.3; He II A6560.0 and Ha 
A6562.8; the [N I] pair at A5197.9, 5200.3; as well as to distinguish nebular He II A4685.7 emission 
from broader stellar emission (when present); and to distinguish nebular from atmospheric emission 
(by virtue of the radial velocity of the targets). Thus in principle it is possible to make use of 
additional emission lines for diagnostics and to determine more accurate ionic abundances from 
lines that are often partially blended. Deeper exposures offer the possibility of employing less 
commonly used diagnostics, such as the faint, temperature-sensitive auroral lines of [O I] A5577.3 
and [N II] A5754.6, or the density-sensitive doublet of [Ar IV] AA4711. 3, 4740.2. Deeper exposures 
also increase the chance of detecting higher ionization stages of important elements, e.g., [Ar IV] and 
[Ne IV], which decreases the need for (or magnitude of) ICFs in computing elemental abundances. 
In practice, though, the low abundances of the SMC PNe mean that the emission lines of most 
ionic species are intrinsically fainter relative to H/3 than their Galactic counterparts, so that deeper 
exposures are needed even for standard abundance determinations. Finally, the larger wavelength 
coverage, particularly the IR, allows us to determine abundances of ionization stages not observable 
in the optical, or to determine ionic abundances with greater accuracy, since they are far less affected 
by interstellar extinction and have little dependence on electron temperature. 



2.1. Optical Spectra 



The observations were obtained with the ESO Multi-Mode Instrument (EMMI) (Dekker, De- 



labre & Dodorico 1986) on the 3.6-m New Technology Telescope (NTT) on La Silla, Chile in 2000 
and 2003. We observed our targets with both long-slit and echelle spectroscopy to obtain the nec- 
essary spectral coverage and resolution with the minimum number of configurations. The long-slit 



spectra were obtained with the blue arm, a sht length of 19'/5 which has a spatial resolution of 
~ (/.'SG pixel"^, and grating #3 which has a spectral resolution of R « 3000. We selected a wave- 
length range of 3500-3900 A for the 2000 observations, but adjusted this to 3700-4100 A for the 
2003 observations to include a small wavelength overlap between the blue and red spectra. The 
echelle spectra were obtained with the red arm using grating ^9 and the cross-dispersing grism 
7^3. The 2003 observations were obtained with an updated red CCD array with 2 chips, each with 
2 independent amplifiers. The space between the detectors produces a gap in the spectral coverage 
at 4900-5000 A. The dispersion is ~ 0.12 A pixel"^ (~ 0.21 A pixel"^ for the spectra taken in 
2000); grating ^^9 gives a spectral resolution oi R ^ 8400 (7700), corresponding to ss 0.6 A near 
the center of the wavelength range. 

The observing log is presented in Table [T} The log includes the target name using the catalog 
designations of [Sanduleak, MacConnell &: Phillipl ( [1978^ and [Morgan fc Good] ( [l992l ), the UT date 



of the observation (s), the exposure times for each arm of the instrument, the width of the entrance 
slit, and the wavelength range covered. Note that most bright MCPNe are smaller than I'.'O, so 
that the slit size could be adjusted for the seeing and still intercept most of the nebular emission 
with little loss of spectral resolution. Because the airmass for our exposures ranged from about 1.2 
to 1.5, we aligned the entrance slit with the parallactic angle for the beginning of the initial pair of 
exposures for each target. Even with this precaution, the parallactic angle changed by as much as 
20° during a visit; regrettably, this telescope/instrument combination does not offer atmospheric 
dispersion correction, nor does it automatically rotate the slit to follow the parallactic angle. 



The long-slit data were calibrated following standard techniques, using IRAIrJ software (Massey 



1997). The bias and flat-field frames were combined by averaging them, with artifact rejection. The 
spectra of the targets were then combined with artifact rejection, bias subtracted, trimmed, and 
fiat-field corrected. The task doslit, within the noao.imred.specred package, was used for spectral 



extraction, and wavelength and flux calibration (Massey et al. 1992). The extracting aperture 
was 3'.'6, which was more than enough to ensure that the entire spatial extent of the nebula was 
included. 



The echelle spectra were calibrated following the steps outlined in Willmarth & Barnes (1994). 
The target spectra, standard star spectra, flat fields, and wavelength calibration files were first bias 
subtracted using an average bias frame, then trimmed to remove the overscan as well as the bluest 
orders that were too faint to trace. The photometric uniformity of the target spectra, standard 
star spectra, and wavelength calibration files were then corrected with the normalized fiat. The 
extraction aperture for the targets was 3'.'5; sky emission was extracted on either side of each order 
with aperture sizes of 3'.'2 (for the 2003 observations) or 2'.'2 (for the 2000 observations), which 
were separated from the target apertures with a gap of ?» 1". Th-Ar comparison arcs provided an 
accurate wavelength calibration. A correction for atmospheric extinction was made using average 

^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



values for the observing site, and a response function was constructed for each order using standard 
star spectra obtained on each night. The flux-calibrated standard star spectra were then checked 
for order-to-order discontinuities, and the extraction and sensitivity function determination were 
refined as necessary to ensure an accurate calibration of the target spectra. 

Although the long-slit and echelle data for a given target were usually obtained consecutively, 
the blue arm (long-slit) spectra needed to be placed on the same intensity scale as the echelle (red- 
arm) spectra. This scaling is necessary because the spectra were taken in conditions that were not 
always photometric or with unstable seeing, so that grey extinction, slit losses, and other factors can 
in principle be different. Our original intent was to use the emission lines common to both spectra 
to establish the scale factor. Unfortunately these lines are weak and few in number, and more 
importantly the sensitivity in the bluest orders of the red arm is low and declining so steeply with 
decreasing wavelength that scale factors so derived were unreliable. Comparisons with published 
data for this purpose proved unreliable as well (see §2.2). In the end we derived the blue/red scale 
factors using the H Balmer decrement, such that the H7, H5, and He intensities were in good 
agreement with that predicted from the /(H/3)/I(Ha) ratio, and the derived electron temperature 
and density. The adopted blue/red scale factors are given in the last column of Table [Tl and are 
usually within 10% of unity. 



2.2. Optical Emission Line Intensities 

We measured the emission line intensities in all the NTT spectra using the IRAF splot routine. 
We directly integrated the flux above the local continuum wherever possible, except for close, 
nearly blended lines where we cross-checked the individual line fluxes with Gaussian deblending. 
The physical size of the red-arm detector is sufficiently large that some emission lines appear on 
opposite ends of two adjacent echelle orders. In these cases we found that the fluxes measured 
in each order agreed to within 10% in most cases; ultimately we selected the flux from the line 
closest to the center of the order, both because of the higher signal-to-noise ratio and because the 
order trace was better defined near the order centers. The level of agreement between orders sets 
a lower limit to the uncertainty of all the fluxes. The flnal nebular emission line intensities are 
presented in Table ^ after application of the red/blue scale factors but uncorrected for interstellar 
reddening. In successive columns we give the identification of the ion responsible for the emission, 
the rest wavelength, and the relative line intensity for each target, normalized to I(H;S)=100. Note 
that the intensity of [O III] A4959 is missing in most cases: this line falls in the gap between the 
detectors in the echellograms obtained in 2003. 

Every one of our targets has been observed spectroscopically by others, although usually with 
some combination of lower resolution, shallower exposures, or lesser wavelength coverage. We have 
compared our optical emission line intensities with those published by several other investigators: 
Monk, Barlow k Clegjl ( |1988[ ), [Meatheringham fc Dopita| ( |1991a|b[ ), [Vassihadis, et al.| ( |1992[ ), 



Leisy k Dennefeld (1996 2006), Idiart, Maciel k Costa (2007, hereafter, IMC07), as well as those 



from the HST /STIS spectra of Stanghellini, et al. (2003). The comparison data to which we 
attach the most weight are the HST /STIS shtless spectra, owing to their superior absolute and 
relative calibration. But even these STIS data have limitations, particularly for targets larger 
than ~ 1'.'2 in size because of the spatial overlap between Hq and nearby lines, and because fluxes 
from extended targets generally had lower signal-to-noise ratios. Therefore the STIS fluxes from 
MG 8 and SMC 11 have been excluded from the following comparison. Figure [T] shows moderately 
good agreement between our fluxes and those in the literature, with the best agreement for strong 
lines and an increase in dispersion at fainter flux levels. A few trends are worth noting, however. 
First, the intensities from Meatheringham & Dopita ( 1991a|b ) are systematically brighter than 
those presented here when the intensity we measure is < 5% of /(II/3), with a dispersion that 
is quite large compared to their quoted errors. Second, the intensities from |Leisy &: Dennefeld 



(1996, 2006) are 17% fainter than ours when the intensity we measure is less than that of H/3, 



although the agreement for brighter lines is excellent. Third, the intensities published by IMC07 
show a large dispersion relative to ours, well beyond their quoted errors, and trend systematically 
brighter by ~8% when the intensity we measure is less than ~5% that of II/3. It should be noted, 
however, that IMC07 adopted mean emission line intensities that were averaged over available 
data in the literature. This complicates the comparison because each emission line was determined 
from a different combination of observations and literature values. Thus, systematic errors, while 
undoubtedly present, are very difficult to characterize. 

The excellent agreement between the line intensities presented here and those from the HST /STIS 
data (within 3% in the mean) suggests that our intensities are the more reliable ground-based mea- 
surements. Based upon the level of internal consistency in our line intensity measurements and 
the above comparison of these intensities with values in the literature, we assign an uncertainty to 
our data of 10% for intensities brighter than 10.0 (on the scale of /(H/3) = 100 prior to correction 
for interstellar reddening), 20% for intensities between 2.0 and 10.0, and 40% for intensities fainter 
than 2.0. Intensities fainter than 0.5 probably have uncertainties approaching 100%. 



2.3. Infrared Emission Line Intensities 



Infrared spectra of several of our SMC PNe were obtained with the Infrared Spectrograph 



(Houck et al. 2004) on-board the Spitzer Space Telescope between 2005 July and November, as 



a part of General Observer program 20443. The spectra from this shallow survey are all low 
dispersion (R~ 100) covering the 5-38 //m range, and of relatively short duration (less than 2 min 



total exposure time per spectral region). Stanghellini, et al. (2007) describe the observing program 
and the data reduction techniques. The nebular emission line fluxes were measured with a Gaussian 
line-fitting routine in the SMARTrl package. Many expected emission lines were either below the 



^ SMART was developed by the IRS Team at Cornell University and is made available through the Spitzer Science 
Center at Caltech. 



detection limit, or were embedded within strong solid-state features. Where possible, we derived 
upper limits by computing the flux of a fictitious Gaussian emission line with a width typical of 
narrow emission features at the expected wavelength, and a height equal to three times the RMS 



noise in the local continuum. Our data were supplemented with emission line fluxes from Bernard- 



Salas et al. (2008), whose similar GTO program included Spitzer /IRS spectra of two targets in 
our sample. The emission line intensities from all the Spitzer spectra, which included all targets 
except for MG 8, MG 13, and SMP 9, are presented in Table [3J In successive columns we give the 
identification of the ion responsible for the emission, the rest wavelength in microns, the reddening 
function /(A), and the relative line intensity for each target, uncorrected for interstellar reddening 
and normalized to /(H/3)=100. For subsequent analysis we correct the intensities for interstellar 
reddening in the standard way by multiplying them by lO'^-' *- \ where c is the logarithmic extinction 



at H/?. The log of the absolute H/3 fiuxes, given in the first row, are from Stanghellini, et al. 



(2003). Although there are no independent IR emission line flux measurements for comparison 
with these data, for the purpose of computing abundance uncertainties we adopt our statistical 
flux uncertainties, which are of order 10% for most lines. 



3. Chemical Abundances 

3.1. Physical Diagnostics 

Determining the interstellar extinction and the physical diagnostics (gas density and tempera- 
ture) is the first step in deriving the elemental abundances of the nebular gas. We have adopted the 



interstellar extinction constants from the accurate spectrophotometry of Stanghellini, et al. (2003), 
with the exception of those for MG 13, for which the Ha flux is not available, and SMP 11, where 
the angular extent of the nebula is large enough that the HST fluxes at Ha and [N II] AA6548, 6583 
suffered considerable overlap. In these cases we derived the extinction constant, c, from our mea- 
surements of the Balmer decrement in the NTT spectra, assuming an intrinsic /(Ha)//(H/3) ratio 
of 2.85 for case B recombination. We used ratios of collisionally excited emission lines and the 



nebular software package (Shaw Sz Dufour 1995; Shaw et al. 1998) to determine the nebular 



electron densities (A^'e) and temperatures (Te). Specifically, we constructed traditional diagnostic 
diagrams for each object, plotting curves in the N^, Tg plane that were consistent with the ob- 
served (dereddened) ratios of [O II] /(A3726.0)]//(A3728.8), [Ar IV] /(A4711.3)]//(A4740.2), and 
[S II] /(A6716.5)]//(A6730.9) for density; and [N II] [/(A6548.0)+/(A6583.4)]//(A5754.6) and [O III] 
/(A5006.8)/I(A4363.2) for temperature. A few other diagnostic curves were plotted as well, includ- 
ing the optical temperature diagnostics [O II] /(A3727)]//(A7325) and [S II] /(A4072)]/I(A6725), but 
in practice the uncertainties in these ratios are too large to be useful for this purpose. They were in- 
stead used as a check on the accuracy of the flux calibration and the extinction correction, for which 
the long wavelength baseline is a strong advantage. In the infrared, it is possible to determine N^ 
from the IR ratios of [S III] /(A18.68 /im)]//(A33.64 /um), and [Ne V] /(A14.29 /um)]/I(A29.23 /im), 
which are sensitive to densities of a few xlO'^ to r^ 10^ cm~^. In practice, the red-ward line 



fluxes are usually upper limits, so at best we can only use them here to exclude the possibility 
of moderate- to high- ionization regions with densities > 10^ cm~^. Temperature indicators from 
ratios of optical to IR emission lines can also be used in principle: those available here are [S III] 
/(A6312.1)]//(A18.68 iim), and [Ne III] /(A3868.7)]//(A15.56 //m). In practice the [S III] A6312.1 
line is too weak for a reliable temperature, although the median of the differences between the Tg 
derived from this ion and that from [O III] is less than 100 K. The [Ne III] lines, while strong, suffer 
from a systematic discrepancy between the IR and the optical emission lines which is discussed 
in §3.3.1. Therefore, the temperatures derived from IR lines are not sufficiently reliable for this 
purpose, but do serve as another check on the relative flux calibration and reddening correction. 

The derived physical diagnostics are presented in Table [4j where the first two columns give 
the target (common) designation and the adopted extinction constant. The next three columns 
give iVe as derived from three diagnostic ratios, with the adopted density in the following column. 
Note that the adopted density is generally taken from the [O II] ratio because the signal-to-noise 
ratio is always higher in these lines; in addition, the [S II] ratios are in some cases near the high 
density limit, meaning that a small uncertainty in the flux ratio yields a very high uncertainty 
in the inferred density. The [Ar IV] lines are extremely weak and are therefore unreliable as a 
density diagnostic. The next two columns give Tg as derived from the available diagnostics; formal 
uncertainties were also determined from the uncertainties in the flux measurements, which were 
taken from the prescription in §2.2. In all but one case the error bar of Te([N II]) overlaps with 
Te{[0 III]), and in two cases Te([N II]) > re([0 III]). In general Te([N II]) would be expected to be 
more representative of the low-ionization zone, which would contain most singly-ionized species. 
The weakness of these lines in many objects, which is largely a result of a low N abundance, 
means that they seldom yielded Te with a sufficiently small uncertainty to be useful. The adopted 
temperatures, given in the last two columns in the Table, are taken from [N II] for the low-ionization 
region (when reliable) , and from [O III] for the medium-ionization region; we adopt Tg ( [O III] ) when 
Te([N II]) is not available or when it exceeds Te{[0 HI])- 



3.2. Ionic Abundances 

The abundance for He"*" was derived from the A5875.7 recombination line. For these calculations 
we used the electron temperature as derived from [O III] , and the case B recombination coefficients 



from Porter et al. (2005) which makes use of the latest radiative and collisional data. The abundance 
from the recombination lines at A4471.5 and A6678.2, was computed as well, but these lines are 
a factor of ~ 3 weaker (and typically only a few percent of H/3), with correspondingly greater 
uncertainty. The average of the abundance as derived from A4471.5 and A6678.2, where available, 
generally agreed well with that from A5875.7. The He"*"^ abundance was derived from the A4685.7 



recombination line, ^£([0 III]), and the recombination coefficients as tabulated in Osterbrock & 



Ferland ( 2006 ) . The results are given in the first two rows of Table 5 for each nebula. 



The ionic abundances for collisionally excited lines were derived using the nebular software 



10 



from the physical diagnostics determined above — where Te for the appropriate ionization zone was 
used (see Shaw k. Dufour 1995 Shaw et al.||1998| ) — and the emission hne intensities from ions of 
N, O, Ne, S, and Ar as compiled in Tables [2] and [3] and corrected for reddening using the Prevot 



et al. (1984) extinction curve for the SMC. We note that the nebular package has been updated 



with more recent energy levels and transition probabilities for various ions (see the compilation by 
Badnell, et al. 2006, for references), including N+^, 0+^, 0+'^, Ne+^, Ne"*"^, and Ne"*"^ which are 
of interest here. The ionic abundances relative to H^ for collisionally excited ions are presented 
in Table [5] for each nebula. The abundances as computed from the mid-infrared lines are listed 
separately and annotated with "IR" in the first column, next to the ion identification. For some 
ions the abundance may be computed from either the optical or IR lines, which provides a useful 
comparison that will be discussed below. The adopted ionization correction factor (see below) is 
also listed for each element following the abundances for the observed ions. 



3.3. Elemental Abundances 



The elemental abundances of He, N, O, Ne, S, and Ar are presented in Table [6j based on the 
analysis presented in this section. The recombination lines of both ions of He are present in the 
optical spectra, so the elemental abundance is simply the sum of the abundances derived from each 
ion. He"*" is detected in all the PNe studied here, but often He"*"^ is not. We note that the emission 



from He II A4685.7 in MG 8 is entirely of stellar origin: it is incorrectly attributed by Vassiliadis, 



et al. ( 1992 ) as nebular emission. Further, Vassiliadis, et al. ( 1992 ) give only an upper limit to the 



intensity of He I A5875.7 of < 5, yet it is clearly present in our spectrum with an intensity of 14.4. 



Thus, the He abundance derived by LD06 based on the Vassiliadis, et al. (1992) data is incorrect 



3.3.1. Ionization Correction Factors 



We derived total elemental abundances for N, O, Ne, S, and Ar, relative to H+, by summing 
the abundances of the observed ions and correcting for unobserved ionization stages using ICFs. We 



generally adopted empirical relations for ICFs from Kingsburgh & Barlow (1994) for all elements 



except S, where we instead adopted the recipe from Kwitter &; Henry (2001). However, some of 



the IR emission lines originate from ionization stages that are not observable in the optical, so that 
ICFs are sometimes not needed or, in the case of upper limits, can be confirmed. The adopted 
ICFs for each element are given in Table [5j and the log elemental abundances are given in Table [6] 
normalized to log H"*" = 12.0. 

The first two ionization stages of oxygen (O"*" and O^^) are detected in the optical spectra for 
all of the PNe in our sample. In the IR spectra O"^'^ is detected in four objects, and tight upper 
limits are established for an additional seven objects. In these cases it is clear that the O abundance 
is properly obtained from the straight sum of the observed ionic abundances (i.e., ICF(O) = 1.0), 
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since higher ionization stages do not contribute significantly. For the three remaining objects where 
no IR spectrum is available, ICF(O) is computed from the He ionic abundances using Eq. A9 from 



Kingsburgh & Barlow (1994) 



For nitrogen, the only collisionally excited emission lines in the optical spectra are of N^ (very 
weak recombination lines of N III are detected in some objects), and no N lines are observable in 



our IR spectra. Although Leisy & Dennefeld (1996) published UV spectra from lUE, the upper 



limits for N"*"^ , N"*"^ and N"*"^ do not set useful constraints on the abundances of these ions in any of 



our targets. Therefore ICF(N) is derived from Eq. A2 of Kingsburgh & Barlow (1994). Note that 
ICF(N) is uncomfortably large for MG 13 and SMP 23. Therefore only a lower limit for N is given 
in Table p^ for MG 13 based on the N"*" abundance alone (although N/0 can still be derived from 
N"*"/©^); for SMP 23 the N abundance is especially uncertain, being estimated from < 1% of the 
total, although no recombination lines of N III were detected. 

Two ionization stages of neon (Ne^^ and Ne^^) are observable in the optical, while four stages 
(Ne^, Ne"*"^, Ne^^, and Ne"*"^) are observable in the IR. Most of these lines are faint, however, and 



in many cases only upper limits are available for stages other than Ne"*"^. Following Kingsburgh & 



Barlow (1994), we assume the Ne"*" contribution is negligible, and set ICF(Ne)=1.5 when Ne+ is 



observed but Ne"*"^ is not; and use their Eq. A28 if only Ne"*"^ is observed. In all cases the ionic 
abundances from the IR lines are used where available, a point we discuss in §4. We also note 
that the upper limits for undetected ions are all consistent with the adopted ICFs. Curiously, the 
abundances of Ne^^ as derived from the IR lines are roughly 50% greater than those derived from 
the optical lines, the latter being derived entirely from [Ne III] A3868.7 emission. We discuss the 



fidelity of the optical measurement in §4.1, but note here that Pottasch & Beintema (1999) found 



a similar effect, though with a much larger discrepancy between the optical and IR. 

As with oxygen, the first two ionization stages of sulfur (S^ and S^^) are observable in the 
optical spectra for all of the PNe in our sample, although S"*"^ is not actually detected in two of the 
objects. In the IR spectra, S"*"^ is detected in six objects, and useful upper limits are established 
for an additional five; S"*"^ is detected in nine objects, and useful upper limits are established for 
two others. A comparison of the S"*"^ abundance from the optical and the IR shows generally 
good agreement, and upper limits in the IR equal or exceed those derived from the optical lines. 
Once again we adopt the S^^ abundance from the IR lines where available. The S'^'^ abundances 
and upper limits shed light on the computed ICFs, which we derived from Eq. 4e of K witter"^ 



Henry (2001) where necessary. For five objects (SMP 8, SMP 11, SMP 13, SMP 24, and SMP 27) 



the abundance from the simple sum of the first three ionization stages (which we adopt, setting 
ICF(S)=1.0) is within ~ 11% of that derived from (S++S+2)xICF(S). For SMP 18 the upper limit 
to S+^ shows that the ICF is within 10% of the correct value. Of the remaining four objects, three 
of them (SMP 14, SMP 17, and SMP 19) have relatively high ionization, with the abundance of S^^ 
exceeding that of S^^ and suggesting that contributions to the total S abundance from even higher 
ionization stages may be significant. In these cases, the S abundance derived from ICF(S) is even 
lower than the S abundance derived from the simple sum of the observed ions; we adopt the simple 
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sum as a lower limit. Finally, S^^ is not detected in the optical spectrum of SMP 23 and only an 
upper limit is available from the IR spectrum, yet S"*"^ exceeds S+ by an order of magnitude. In 
this case we adopt the straight sum of the abundances from the observed ions as a lower limit. 

All of the relevant ionization stages of Ar are observable in the optical, except for Ar^ which 
is only observable in the IR. The Spitzer /IRS spectrograms were too shallow to detect the [Ar II] 
6.99 una line, so only upper limits to Ar"*" could be computed. In all cases these limits are much too 
large to be useful for constraining the derived ICF(Ar). Therefore the prescription for ICF(Ar) is 



taken from Eq. A30 of Kingsburgh Sz Barlow (1994) to correct for its absence. Ar"*" is also observed 



in the IR, but only upper limits are available for these targets; the limits are all consistent with 
the ionic abundances determined from the optical spectra. 



3.3.2. Uncertainties 

We computed the uncertainties for each element abundance by determining the effect of the 
uncertainty in Tg, and adding in quadrature the effect of uncertainties in the emission line intensities 
that were quantified in §2. The uncertainty for abundances determined solely from optical emission 
lines is ~ 0.10 dex, owing to the strong dependence of abundance on Tg, while for the IR lines 
it is ~ 0.05 dex. Thus, the magnitude of the uncertainties in the elemental abundances varies 
considerably, depending upon the wavelengths of the transitions for the dominant ions. Broadly 
speaking the uncertainties are 0.1 dex or better, but may be higher in individual cases as noted 
above. No allowance has been made for the statistical or (the likely larger) systematic uncertainty 
in the derived ICFs, but since the corrections are generally of modest size (except for N) these 
uncertainties would not dominate the errors. Generally speaking, the magnitude of our adopted 
ICFs is comparable to those in the literature (see, e.g., Kingsburgh &: Barlow|1994 Kwitter Sz Henry 



2001, Henry, Kwitter Sz Balick 2004) in those cases where only optical emission lines are available 
for analysis. The main exception is for MG 8, which will be discussed separately below. The 
relatively large ICF(N) for most of our objects results from the availability of only one observable 
ionization stage (N"*") in the optical and IR, the moderate- to high-ionization of most of this sample, 
and from the relatively low enrichment of N in these non-Type I PNe. 



3.4. Comparison to Published Abundances 

Our sample of objects presents an unusual opportunity to examine in detail the sources of 
uncertainty in the determination of elemental abundances using direct methods, and to explore the 
discrepancies between different studies. All of the objects in this sample have been observed by at 
least two, and usually a few observers. They are all angularly small, so the entire nebula was most 
likely included in the spectrograph aperture, and therefore ionization stratification and differential 
effects of dust are not factors in the comparison. All of these objects have been observed with at least 
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one space-based spectrograph: most have been observed with both HST/STIS and Spitzer /IRS. 
Thus, the absolute fluxes of key emission hues, and the relative fluxes of many other lines over a 
broad wavelength range, are well determined. It is especially useful to compare abundances that 
are derived from optical lines to those derived from the IR, both as a check of consistency and as 
a means of evaluating the veracity of the ICFs. The advantages of deriving abundances from the 



IR flne structure lines have been enumerated before (see, e.g., Bernard-Salas et al. 2001). Briefly, 
the IR emission lines are stronger than the collisionally excited emission lines in the optical, they 
are not greatly affected by interstellar extinction, the abundances derived from these lines are little 
affected by uncertainties in Te, and more ionization stages can be observed than are available in the 
optical alone. To this we would add that the IR emission line strengths measured from Spitzer /IRS 
spectra are more accurately calibrated than those from most ground-based spectra (optical or IR). 
In §3.3.1 we compared the elemental abundances as determined from optical emission lines and 
ICFs to those derived from a combination of optical and IR lines, without ICFs. This idea is not 



new and appears to have been flrst suggested by Pottasch & Beintema (1999), but here we are able 



to make comparisons for a greater number of ICF prescriptions. 

We compared our elemental abundances in Table [6] to those presented in the large compilations 
of LD06 and IMC07 for the same objects. The last four rows of the Table give the median abundance 
(excluding limits) for targets in this sample, the median abundances from LD06 and from IMC07 
for targets in common (excluding highly uncertain abundances), and the mean of H II regions in the 
SMC from Dennefeld (1989). Figure^ shows histograms of the differences in A(X), the logarithm 



of the abundance for element X, over all elements. The resulting distributions are quite broad: 
the half-width exceeds 0.15 dex and large wings extend to 1.5 dex. It is also useful to compare 
the differences in the median abundance for each element between the two studies, which we give 
in the last few rows of Table [6J At the outset we expected an offset of up to +0.07 dex and 
—0.03 dex (for LD06 and IMC abundances, respectively) derived from faint lines because of the 
flux discrepancy identifled in §2.2. In fact, we find our abundances for He are higher by 0.16 dex. 
While the origin of this discrepancy is not entirely clear, part of the reason may be our use of the 
more recent calculation of the He I recombination coefficients in the Case B limit by |Porter et al 



(2005). For N, our median abundance is lower than those of LD06 by 0.34 dex, although if the 
values noted by LD06 as poorly determined are excluded, then our median abundance exceeds that 
of LD06 by only 0.14 dex for PNe in common, which is within their quoted uncertainty. Our median 
abundance is less than than that of IMC07 by 0.09 dex, which is within the quoted uncertainties. 
We believe our N abundances are the more accurate, owing mostly to improved accuracy of the 
emission line intensities. Our median O abundance is lower by 0.07 dex and 0.06 dex than that of 
LD06 and IMC07, respectively, which is within the quoted uncertainties in both studies. However, 
this similarity belies the much larger range of values, which will be discussed in §4.2. 

Our median Ne abundance is higher by 0.22 dex than that of LD06, but this is almost entirely 
due to our use of abundances from the IR lines; had we used only the optical lines, the difference 
would have been near +0.04 dex. Our median Ne abundance agrees well with that of IMC07, to 
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within 0.05 dex. It is not clear why the abundances from the optical and IR emission lines are so 
discrepant for this element. Two potential sources of error, the flux calibration and the extinction 
constant, were validated by the consistency of the ratios of [O II] /(A3727)//(A7325) and [S II] 
/(A4072)//(A6723) with the adopted extinction constants and Tg, and by the consistency of the 
ionic abundances as derived from the optical and IR for other elements. Since a similar discrepancy 



was noted by Pottasch Sz Beintema (1999), we suspect the problem may lie with the supporting 
atomic data. Our median S abundance is 0.66 dex and 0.61 dex lower than LD06 and IMC07, 
respectively. Although this difference is exceptionally large, in the first case it is based on only 
the 3 objects in common that were not excluded from the LD06 study as unreliable. Furthermore, 
the [S II] line fluxes are weak: the strongest in this sample are ~ 10% of H/3, and the weakest are 
~ 1%. Thus many of the S abundances are subject to the systematic errors of the type described 
in §2.2. We are confldent in our S abundances because the combined optical and IR spectra cover 
most ionization stages of S, and that ICF(S) as computed from the optical emission lines is in 
excellent agreement with that inferred from the abundances as computed from both optical and 
IR lines. Finally, our median Ar abundance is only 0.08 dex and 0.06 dex lower than that of LD06 
and IMC07, respectively, and well within the quoted uncertainties for this element. 



4. Discussion 

4.1. Errors in Elemental Abundances 

The uncertainty in the determination of the elemental abundances using direct methods de- 
pends upon a variety of primary and secondary factors, including the uncertainty in the observed 
emission line intensities, uncertainties in the extinction, Tg and Ng, and in the determination of the 
ICFs, which themselves depend upon abundances of other ionic species. It has often been argued 
that the majority of abundance uncertainties are predominantly statistical in nature, and are gov- 
erned by shot noise in the emission line intensities. While it is true that some ionic abundances are 
derived from very weak emission lines, the majority are derived from lines with intensities that are 
typically greater than 10% that of H/3, with corresponding quoted uncertainties of order ~ 10%. 
Yet Figure p^ shows that the abundance discrepancies found here are ~ 0.15 dex (40%) or more. 
We argue that the magnitude of the discrepancies identifled here includes a significant component 
of systematic errors in emission line intensities, which when propagated through the analysis result 
in greatly magnified errors in elemental abundance determinations. 

Systematic errors in emission line intensities can have a variety of causes, including flawed 
observing technique and poor flux calibration. If, for example, a ground-based spectrograph is 
configured to obtain a large wavelength coverage with a narrow slit, the resulting observations of 
compact targets may be compromised by color-dependent slit losses due to atmospheric dispersion 



even at modest airmass (Filippenko 1982) unless care is taken to mitigate the effect during the 



observation. The effect of such slit losses, which are often most severe at the blue end of the spectral 
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range, can be manifested in inferred interstellar extinction constants that are systematically too 
large. Figure p^ compares the extinction constants (expressed as c, the logarithmic extinction at 



H/3) used here with those measured by Monk, Barlow Sz Clegg| (1988), Meatheringham Sz Dopita 
( 1991a|b ), Leisy & Dennefeld (1996, 2006), and IMCOTJ^ Errors of this sort, or in flux calibration 
can be detected using sanity checks with emission lines that span a large range of wavelength, such 
as ensuring that the [O II] /(A3727)]//(A7325) and [S II] /(A4072)]//(A6725) ratios are consistent 
with the adopted c and Tg. Generally, errors in c smaller than 0.2 do not cause large errors in ionic 
abundances derived entirely within the optical band, but errors this large may result in substantial 
discrepancies between ionic abundances derived from UV and IR emission lines. Errors in the 
extinction constant are unlikely to cause significant errors in the determination of T^, as the most 
commonly used ratios (from O III I(A5006.8)//(A4363.2) and [N II] [/(I(A6583.4)]/I(A5754.6) span 
relatively modest wavelength baselines. 

The ionic abundances depend linearly on the observed emission line intensities. However, the 
dependency on Te is either minor (for transitions in the IR) or exponential (for optical and UV 
transitions) . It is this exponential dependence of the optical emission lines on T^ that can inflate 
modest systematic errors in fluxes to large errors in ionic abundances. Systematic errors of only 
10% in emission line intensities can lead to errors in Tg of ~ 1000 K, depending on the diagnostic 
used. For temperatures in the range considered here, an error of 1000 K can translate to an error 
in ionic abundances as determined from optical emission lines of ~ 0.08-0.15 dex. Judging by the 
magnitude of the discrepancies between Te as determined here and those of LD06 and IMCOTq 
shown in Figure [4j it is not surprising that the elemental abundance discrepancies can be much 
larger than the quoted statistical errors for the individual line fluxes. Other factors can inflate 
abundance discrepancies as well, such as errors in the supporting atomic data and the adopted ICF 
recipes. The former can lead to discrepancies of 30% if the latest atomic data are not used, while 
in practice there is considerable agreement on the ICF recipes. 



4.2. Abundances as Tracers of Prior Stellar Evolution 



The small number of objects in this sample of PNe makes it difficult to gain deep insight into 
nucleosynthesis and dredge-up processes of AGB stars in the SMC, but a few important conclusions 



can be drawn. Villaver, et al. (2004) determined the central star masses of 10 of the PNe in this 
sample, and with two exceptions they are ~ 0.59 Mq, implying moderate- to low-mass progenitors. 



In the end, IMC07 adopted the mean of extinction constants from several sources. 

*The largest apparent discrepancy in Te was for MG 8, where IMC07 derived 23,900 K compared to our 13,000 K. 
However, they have adopted an intensity of 5% of H/3 for the auroral emission line of [O III] A4363, which appears to 
have originated from LD06. This emission line was in fact listed by LD06 as an upper limit to the intensity from the 



spectrum published by Vassiliadis, et al. (1992 1. The ratio of this intensity to the intensities of the nebular [O III] 
lines from Vassiliadis, et al. ( 1992| would yield the erroneous Te that was used by MC07. Since this appears to be 
an error, rather than a difference in measurement or technique, we do not include that high Te in Fig. [4| 
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This inference is reinforced by the PN abundances: Figure [5] shows N/0 vs. He for our sample of 
SMC PNe, along with the data from LD06 and IMC07 for the same objects. This diagram has 
traditionally been used as a diagnostic to discriminate between Type I and non-Type I MCPNe 



( Torres- Peimbert Sz Peimbert 1997). Based on morphological and kinematic data, Type I PNe in 



the Galaxy are associated with bipolar shapes and higher mass progenitors. Chemically, they show 
substantial enrichments of N and lesser enrichments of C relative to O, which is consistent with 
theoretical predictions for hot-bottom burning (HBB) in the most massive PN progenitors at the 
end of AGB evolution (see, e.g., Marigo, et al. 2003| ). Non-Type I PNe show lesser enrichments 
in N/0 and He, substantial enrichments of C, and typically have round or elliptical morphologies. 
There are clearly no Type I PNe in our sample, although all of the nebulae except SMP 11 show 
high N/0 and high He and N abundance relative to H II regions. This is consistent both with 
the nebular morphology (again, except for SMP 11), and with the high abundance of C in the 9 



objects where it has been measured (Leisy Sz Dennefeld 1996: Stanghellini, et al. 2009). MG 8 is 



an especially interesting object in this context in that Villaver, et al. (2004) derived a central star 
mass of 0.88 Mq, implying a progenitor Main Sequence mass of at least 3.5 Mq according to the 



models of Vassiliadis & Wood (1994). This is close to the minimum mass that is thought (at SMC 



metallicities) to be required for hot-bottom burning (HBB) during the prior AGB evolution, which 
would have converted C to N very efficiently. The comparatively modest N abundance in MG 8 
sets a useful constraint on the maximum progenitor mass in the SMC that still does not undergo 
HBB during AGB evolution. 

Further insight into nucleosynthetic processing in the progenitor AGB stars can be gained from 
Figure [6l which shows N/0 vs. O in the nebular gas, and compares our data to that from LD06 
and IMC07 for this sample. The rather large scatter in the LD06 and IMC07 data for both N/0 
and O is greatly reduced in our data, which again illustrates rather dramatically the benefit of 
improved accuracy in abundance determinations. If SMP 11 is excluded, the points are consistent 
with no decrease in N/0 with increasing O. Although our sample is rather small, stellar evolution 



models (Marigo, et al. 2003) predict that the 0-N cycle will not be activated for AGB stars that 



do not experience HBB, which is consistent with the lack of a significant trend (i.e., zero slope) in 
Figure [6] for the progenitor AGB stars in this mass range. 



Planetary nebulae can also probe stellar evolution theory for massive stars (see, e.g., Henry, 
Kwitter Sz Balick||2004 ) by studying the abundances of elements such as Ne, S, and Ar that are not 
predicted to be altered during the course of prior AGB evolution; these abundances would therefore 
reflect that of the interstellar medium when and where the progenitor star was formed. Very often, 
these elements are compared to the abundance of O which is also assumed to be unaltered for stars 
in this mass range. Figure [7] shows the abundance of Ne vs. O, which is well fit with the relation 
log (Ne) = -1.96 + 1.16 log (O), with an RMS of 0.11 when MG 8 is excluded. This fit is in 
nearly perfect agreement with that of Henry, Kwitter &; Balick (2004). The relation also agrees 



with the results of Stanghellini, et al. (2006) for Galactic PNe, also shown in Fig. ItI in the sense 
that their trend line would intercept the current sample of SMC PNe which have a lower average 
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O abundance of roughly half a dex. The tight correlation is strong evidence that Ne and O vary 



in lock-step, as argued by several investigators (Henry, Kwitter Sz Balick 2004 Richer & McCall 



2006 Maciel, Costa & Idiart 2009), and that therefore O is basically unaltered during the course 



of evolution of stars that produce PNe. 

Concerning S, we note that log(S/0)= —1.98 (in the median) for this sample, which is lower 
by ~ 0.3 dex compared to that for SMC H II regions (see Figure^. This is consistent with the 



"S anomaly" identified by Henry, Kwitter & Balick (2004) in Galactic PNe, who suggested that it 



most likely results from erroneous S abundance determinations owing to an inadequate accounting 
of unobserved S^^ in the ICF(S). But since our S abundances are derived mostly from mid-IR 
emission lines, S+^ is detected if it is present in significant quantity. Thus, the S anomaly appears 



to be genuine, and persists even in a low metallicity environment. Henry, Kwitter &; Balick (2004) 
also suggested a similar deficit for Ar, and we find log(Ar/0)= —2.47 (in the median) for our 
sample. This value is also lower by ~ 0.3 dex compared to that for SMC H II regions. The case for 
an Ar deficit is a bit less convincing (see Figure ^ since we rarely detected the very faint Ar^'^ and 
^j.+4 gniiggion lines in the optical spectra. One possible explanation may be that the abundances 
of S and Ar in the interstellar medium of the SMC have changed relative to O and Ne over the 
lifetimes of the progenitors of these PNe. However, we would caution that the source of the H II 
region abundances ( Dennefeldl|1989 ) is some twenty years old, and that a study of a larger sample 
with greater depth, better and more complete wavelength coverage, and updated atomic data is 
perhaps overdue. 



4.3. Unusual Planetary Nebulae 

MG 8 is a very low excitation nebula, but its slightly elliptical morphology leaves little doubt 



of its classification as a PN (Shaw, et al. 2006), although it could be a bipolar viewed pole-on. Its 



spectrum shows broad stellar emission from He II and C IV, indicating the presence of a wind. It 
is unusual in that the abundance of Ne is low by about an order of magnitude for its O abundance 
(see Fig. JTI), yet the Ar abundance is high by about 0.5 dex (see Fig.^ compared to the other PNe 
in this sample. We note, however, that the uncertainty in the O and Ne abundances is relatively 
high for this object, at 0.19 and 0.18 dex, respectively. The ICF(Ne) is higher for MG 8 than the 
other PNe, but the absence of nebular He"*"^ and the modest level of 0"*'^/0 suggests that an even 



higher ICF(Ne) is not warranted. Villaver, et al.| (2004) determined for this object a central star 
mass of 0.88 Mq, Teff = 66,500 K, and log L/L© = 4.33. Although the high mass raises the 
possibility of a higher efficiency for depletion of Ne via ^^Ne(a, n)^^Mg processing in the progenitor 
AGB star, it is difficult to see how it could result in an order of magnitude deficit. The abundances 
of the other elements in MG 8 are typical for the PNe in this sample, although no C abundance is 
available. 



SMP 11 was classified by Stanghellini, et al. (2003) as bipolar, with a significantly larger 



extinction than most SMC PNe. However, its chemical composition is quite typical of the average 
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of H II regions in the SMC (Dennefeld 1989). More specifically, N is not elevated as it is for the 
other PNe, and certainly not as much as one would expect for a bipolar morphology. Given that its 
morphology and excitation is not inconsistent with a compact H II region, its earlier classification 
as a PN is in some doubt. 



5. Conclusions and Future Work 

We have presented new emission line fluxes for several SMC PNe based on deep, high resolu- 
tion optical spectra and low-resolution mid-IR spectra. We have used a variety of techniques to 
assure the accurate calibration of our measurements, and have endeavored to determine meaning- 
ful statistical uncertainties in the derived elemental abundances. Our optical fluxes show excellent 
agreement with extant space-based fluxes, but the agreement with prior ground-based studies shows 
significant, systematic discrepancies especially at faint intensity levels. Such discrepancies propa- 
gate to differences in quantities derived from ratios of emission line fluxes, such as the extinction 
constant, T^, ionic abundances, and ICFs, the cumulative effect of which can result in comparatively 
large discrepancies in derived elemental abundances among authors. 

We derived ionic abundances using IR fine-structure emission lines where possible, which have 
the advantages of high emissivity relative to optical lines, and insensitivity to uncertainties in 
the derived Tg and extinction. Combined with the superior calibration of the space-based IR data 
compared to ground-based data, the derived ionic abundances have substantially lower uncertainty. 
The IR emission lines are important for deriving ionic abundances for several ions that are not 
observable in the optical. They also provide an important check of the ionic abundances derived 
from the optical data for ions in common, where the agreement is very good except for Ne"*"^. 
We were able to validate elemental abundances derived from the combination of optical emission 
lines and traditional ICFs with those derived from a combination of optical and IR lines (or limits 
thereto). We found excellent agreement for O, S, and Ar; the poor agreement for Ne is due to a 
consistent offset in the Ne"*"^ abundance, which may point to a problem with the supporting atomic 
data. In the end, our data have significantly improved the accuracy of the elemental abundances 
for Ne and S in the target nebulae, with considerable improvement in the abundances of O and Ar 
as well, as compared to those in the literature. 

We have explored evidence of nucleosynthetic processes in the progenitor AGB stars, and 
we confirm the expectations from stellar evolution theory that the low-mass PN central stars 
in our sample show no signs that they underwent HBB on the AGB, nor is there evidence of 
destruction of O through the 0-N cycle. One object in this sample, MG 8, was evidently born of an 
intermediate mass progenitor shows peculiar abundances of Ne and Ar. Given its lack of significant 
N enrichment and inferred progenitor mass of 3.5 M©, this object sets a useful constraint on the 
maximum progenitor mass at this metallicity that does not undergo HBB during AGB evolution. 
We confirmed that the Ne and O abundances vary in lock-step, join smoothly with the relationship 
seen in Galactic PNe, and likely are not affected by prior AGB evolution. We also showed that 
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the S anomaly (i.e., low S/0 relative to local H II regions) seen in Galactic PNe appears to be 
genuine, and suggest that such an anomaly might also be present for Ar at least for our SMC 
sample. Finally, we noted that the abundances of SMP 11 closely match those of SMC H II regions, 
which suggests it may not be a FN after all. 

It is important to note the much larger dispersion in the distributions of He, O, N/0, and to 
a lesser extent He, in the LD06 and IMC07 data for objects in common, shown in Figs. [5} [6j and 
[9j The reduction in the scatter in these diagrams is remarkable, and changes the interpretations 
that have been offered in the literature. Certainly our smaller range in O abundance is much 
more consistent with the average for H II regions, and would not require some mechanism for O 
production during post-AGB nucleosynthesis, as LD06 suggested. The lower scatter in Fig.[5]leaves 
no doubt that there are no Type I FNe in our sample, and if the scatter in the N/0 vs. He/H 



diagrams of LD06 and Maciel, Costa Sz Idiart (2009) is dominated by large abundance errors, it 



would imply that the reported lack of segregation into Type I and non-Type I nebulae could largely 
be a result of large errors masking the underlying relationships. Similarly the evidence in Fig.[6]that 
O is not destroyed in AGB stars that produce non-Type I PNe would not have been convincing had 
the LD06 and IMC07 abundances been used for this study. Finally, our lower and more accurate 



S abundances compared to LD06 and Maciel, Costa Sz Idiart (2009) make a more convincing case 
for a genuine S anomaly. 

It is clear that careful attention to the calibration of deep optical spectra, along with near- 
IR spectra make it possible to achieve much higher levels of accuracy in the derived elemental 
abundances of PNe. Higher accuracy is very important for comparing abundances in PNe to the 
predicted chemical yields for the prior phase of AGB evolution. In a subsequent paper we will 
present similar optical and IR spectroscopy for a larger number of LMC PNe, and we will then be 
able to address nucleosynthetic yields from AGB stars more directly, and contrast them with these 
results for the SMC sample. 
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Fig. 1. — Comparison of published emission line intensities (Ref) in common with those reported 



in this paper, with the 1:1 relation indicated {dashed line). Symbol key: Meatheringham & Dopita 



( 1991a|b ) [circles)] Leisy & Dennefeld (1996, 2006[ ) (gg^areg); [Vassiliadis, et al" (1992) (asterisks); 
Monk, Barlow k Clegjl ( |1988| ) ( ^rmng^eg) ; [Idiart, Maciel fc Costa] ( |2007[ ) {inverted triangles); HST 
spectra from Stanghellini, et al. ( 2003| ) {filled diamonds). [See the electronic edition of the Journal 
for a color version of this figure.] 
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Fig. 2. — Discrepancies in elemental abundances between the results of this paper and those of 
LD06 {open boxes) and IMC07 {filled boxes). Plot shows histograms of the differences in the log 
abundance over all elements considered here. 
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Fig. 3. — Adopted extinction constant from this paper vs. those of Monk, Barlow & Clegg (1988), 



Meatheringham fc Dopita| ( |1991a|bD , |Leisy fc Dennefeld|p996l [20061 ), [idlart, Maciel fc Costa] ( [20071 ) , 
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Fig. 4. — Adopted Tg (in 1000 K) from this paper vs. the difference in Tg between those of other 
authors (LD06; and IMC07) and this paper. Symbols as in Fig. [l} Differences exceeding 1000 K 
(i.e., symbols lying outside the grey box) can result in differences in inferred abundance exceeding 
0.10 dex. 
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Journal for a color version of this figure.] 
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Fig. 8. — S/H vs. 0/H for this sample of SMC PNe. Symbols as in Fig. pi 
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Table 6. Elemental Abundances 



log He 



loeN 



logo 



log Ne 



logs 



Ar log N/0 



MG 8 
MG 13 
SMP 8 
SMP 9 
SMP 11 
SMP 13 
SMP 14 
SMP 17 
SMP 18 
SMP 19 
SMP 20 
SMP 23 
SMP 24 
SMP 27 



11.08 ±0.04 
10.99 ±0.04 
11.11 ±0.04 

11.05 ±0.04 
11.01 ±0.04 
11.11 ±0.04 

11.13 ±0.04 

11.14 ±0.04 

11.06 ±0.04 

11.09 ±0.04 
11.14 ±0.04 
11.13 ±0.04 

11.13 ±0.04 

11.14 ±0.04 



7.15 ±0.13 
> 5.12 
7.05 ±0.11 
7.25 ±0.15 
6.52 ±0.09 
7.30 ±0.11 
7.36 ±0.09 
7.38 ±0.09 
7.11 ±0.10 
7.28 ±0.16 
6.95 ±0.12 
7.14 ±0.16 
7.17 ±0.09 
7.17 ±0.12 



8.15 ±0.19 
8.00 ±0.11 
7.88 ±0.11 
8.32 ±0.16 
8.02 ±0.12 
8.06 ±0.11 
8.29 ±0.09 
8.21 ±0.11 
7.90 ±0.12 
8.14 ±0.10 
7.74 ±0.11 
7.93 ±0.11 
8.06 ±0.11 
8.00 ±0.11 



6.39 ±0.18 
7.12 ±0.12 
7.37 ±0.05 
7.51 ±0.12 

6.90 ±0.05 

7.35 ±0.05 
7.65 ±0.05 
7.67 ±0.05 
7.57 ±0.05 
7.62 ±0.05 

6.91 ±0.05 
7.37 ±0.05 

7.36 ±0.05 
7.20 ±0.05 



6.09 ±0.14 

6.11 ±0.05 
6.34 ±0.14 
6.28 ±0.05 
5.96 ±0.16 

> 6.27 

> 6.15 
6.18 ±0.05 

> 6.05 
5.67 ±0.18 

> 5.93 
6.11 ±0.05 
5.84 ±0.05 



6.30 ±0.10 
5.61 ±0.11 
5.43 ±0.12 
5.88 ±0.12 
5.97 ±0.08 
5.46 ±0.11 
5.82 ±0.11 
5.56 ±0.11 
5.67 ±0.10 
5.54 ±0.11 
5.03 ±0.11 
5.41 ±0.11 
5.58 ±0.11 
5.38 ±0.11 



Note. — Abundances of element X expressed as 12+log N(X)/N(H''') 
''Based on N+/0+. 



-1.00 

-0.51^ 

-0.83 

-1.08 

-1.50 

-0.76 

-0.93 

-0.83 

-0.79 

-0.86 

-0.79 

-0.80 

-0.89 

-0.83 



Median 


11.11 


7.17 


8.04 


7.40 


6.11 


5.57 


-0.83 


Median (LD06) 


10.95 


7.03 


8.11 


7.18 


6.77 


5.65 


-1.05 


Median (IMC07) 


10.95 


7.26 


8.10 


7.35 


6.72 


5.63 


-0.90 


Mean H ll'^ 


10.90 


6.46 


7.97 


7.22 


6.32 


5.78 


-1.51 



''Data from the compilation of Dennefeld (19891. 



